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Stationary states in KPZ type growth have interesting short distance prop- 
erties. We find that typically they are skewed and lack particle-hole symmetry. 
E.g., hill-tops are typically flatter than valley bottoms, and all odd moments of the 
height distribution function are non-zero. Stationary state skewness can be turned 
on and off in the 1+1 dimensional RSOS model. We construct the exact stationary 
state for its master equation in a 4 dimensional parameter space. In this state steps 
are completely uncorrelated. Familiar models such as the Kim-Kosterlitz model lie 
outside this space, and their stationary states are skewed. We demonstrate using 
finite size scaling that the skewness diverges with systems size, but such that the 
skewness operator is irrelevant in 1+1 dimensions, with an exponent y s k — —1, and 
that the KPZ fixed point lies at zero- skewness. 
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1. INTRODUCTION 



Crystal surfaces display interesting scaling properties during growth. One of 
the dynamic processes that has been at the center of attention is so-called KPZ 
type growth, named after the Langevin equation 

^ = v + uV 2 h + ^X(Vh) 2 + v (1.1a) 

with uncorrected Gaussian noise 

(77(1-!, txMra, t 2 )) = 2Do( ri - r 2 )o(t! - t 2 ) (1.16) 

studied by Kardar, Parisi, and Zhang [1]. The growth rate v is modified by the local 
curvature of the surface (the f-term), its local slope (the A-term), and random fluc- 
tuations (the r/-term). It is well established by now that many microscopic growth 
processes belong to the KPZ universality class [2-5]. Such microscopic models have 
been investigated by numerous Monte Carlo (MC) simulations and finite size scaling 
(FSS) studies using exact diagonalization of master equations. Moreover, at least 
one model is exactly soluble in 1+1 dimensions (D=l) [6-9]. Intriguing relations 
with different aspects of physics have been established. Eq.(l.l) is equivalent to 
the Burgers equation for randomly stirred fluids (V/i represents the fluid velocity) 
[10]. KPZ growth maps onto to the directed polymer problem [11,5], and relates to 
equilibrium liquid crystal phases [12]. In ID it describes persistent currents in metal 
rings in the context of asymmetric exclusion models [13-17], and is also equivalent 
to the equilibrium statistical mechanics at facet-ridge end-points in 2D crystals [13]. 

The scaling properties of KPZ type growth have been established, but our un- 
derstanding is not yet at the level we would like. It is useful to make a comparison 
with equilibrium critical phenomena. Master equations for microscopic models, like 
the body-centered solid-on-solid (BCSOS) model and the restricted solid-on-solid 
(RSOS) model, play a similar role as Ising type microscopic models in equilibrium 
critical phenomena. KPZ type Langevin equations play a similar role as <p 4 type 
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field theories. There are roughly three levels at which a particular type of scaling 
invar iance can be understood. The first level is to establish empirically (experi- 
mentally and numerically) the existence of scale invariance and universality. For 
equilibrium critical phenomena this was achieved in the early seventies, for KPZ 
type growth only a few years ago. 

The second level is to obtain analytical confirmation of the empirical scaling 
properties. Exactly soluble models confirm the empirical scaling properties of 2D 
equilibrium critical phenomena. Mean field approximations and Landau- Ginzburg 
theory yield the existence of an upper critical dimension and makes it possible to 
carry out controlled renormalization transformations that demonstrate the existence 
of the Ising fixed-point in (p 4 theory. For KPZ type growth, we only reached this 
level in ID. The BCSOS growth model is exactly soluble in ID. The KPZ fixed point 
in general D has eluded us thus far. The existence of an upper critical dimension is 
yet unclear. The e-expansion analysis of eq.(l.l) [1,10,18] describes the reversal of 
the stability of Edwards- Wilkinson (EW) growth (the point A = 0) with respect to 
A at D = 2. It does not yield a fixed point for the KPZ universality class, except a 
strong coupling one in ID. 

The third level is to express the scaling properties in terms of a free field 
theory. In general this is not possible at all. However, for 2D equilibrium critical 
phenomena such a description emerged during the last 15 years. Coulomb gas 
representations [19,20] and conformal field theory [21] provide a full free field theory 
description of virtually all 2D equilibrium phase transitions. There is no assurance 
that this can be generalized to dynamic processes in ID, but we have a good chance 
since the time evolution operators of master equations in 1+1 dimensions resemble 
closely transfer matrices of 2D equilibrium critical phenomena. Generalizations of 
conformal invariance are being considered [22,23], but it is too early to tell whether 
this will work. 

In this paper we present a master equation study of a generalized RSOS growth 
model. In this model, nearest neighbour columns of particles can only differ by 
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dh = 0, ±1. We consider the most general growth rule involving only nearest 
neighbour step configurations. This gives rise to a 5 dimensional phase diagram. 
We address four issues: (a) we describe the global structure of this phase diagram: 
(b) we derive the exact analytic form of the stationary state in a 4 dimensional 
subspace; (c) we point out that the stationary state is generically skewed; and (d) we 
investigate whether stationary state skewness scales in accordance with conventional 
renormalization theory and the concept of universality. 

Somehow, stationary state skewness seems to have escaped everyone's attention 
thus far. The only discussions of skewness in the literature of which we are aware, 
concern temporal skewness [24-26,5]. One possible reason is that almost all previous 
studies involve MC type simulations. In master equation studies, like ours, the 
stationary state and its properties are directly accessible (as the eigenvector of the 
largest eigenvalue of the time evolution operator). Consider the moments of the 
height distribution h x (t) at a certain moment in time, in a finite ID system with 
periodic boundary conditions, h x+ L(t) = h x (t), 



The first moment vanishes since it defines the average surface height h av (t) . The 
second moment is the conventional measure for the width of the interface. The third 
moment characterizes the skewness. In EW type growth all odd moments vanish 
because of particle-hole symmetry at A = 0. KPZ type growth lacks that symmetry. 
On a local level, skewness means that the curvature at a typical valley-bottom is 
larger in magnitude than at a typical hill-top (or the other way around). Fig.l 
illustrates this. It shows the deterministic time evolution of a sinusoid initial state, 
according to eq.(l.l) without noise for positive A. 

It is important to distinguish between stationary state skewness and temporal 
skewness. The skewness in Fig.l for the deterministic KPZ equation is a transient 
phenomenon. This initial state decays to zero, to a stationary state with no skewness 
whatsoever. Noise pumps sinusoid waves at random frequencies into the surface 




(1.2) 
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(from the spatial Fourier transformed perspective) at all time scales. Each mode 
decays with the same sign for skewness. So in the presence of random noise the 
stationary state tends to be rough and skewed. This picture is too simplistic. It 
ignores the non-linearly coupling between the modes. Moreover, the stationary state 
of eq.(l.l) is known exactly in ID. It is Gaussian. All moments W n for n > 2 vanish 
[5,10,27]. 

The BCSOS model is one of the simplest microscopic KPZ type growth rules 
[28], and it is exactly soluble in ID [6-9]. Its stationary state is trivial as well. It 
is the completely disordered state, without any skewness (nor any other structure), 
not only in the thermodynamic limit but also for all finite L. The origin of this 
is a special symmetry of the time evolution operator. We review this in section 2. 
Throughout this paper we will compare our results for the RSOS model with the 
properties of the BCOS model. 

Our study of the RSOS model illustrates that the simplicity of the stationary 
states of eq.(l.l) and the BCSOS model is accidental. Generically, stationary states 
are skewed. Skewness can be turned on and off in the RSOS model. In section 3, 
we derive the exact form of its stationary state in a 4 dimensional subspace. This 
state is simple. Skewness and all other correlations between steps in the surface are 
absent. Steps are placed at random, with not even nearest neighbour correlations 
between them. The only parameter is the step density. Familiar special points, 
such as the Kim-Kosterlitz (KK) model [29], and speical lines [30] lie outside this 
subspace. Their stationary states are non trivial. Stationary state skewness is one 
aspect of this. 

In section 4, we present a mean-field type derivation of the deterministic part of 
the KPZ equation for the BCSOS and RSOS model. This derivation provides more 
insight into the general structure of the RSOS model phase diagram. In particular, 
it identifies a second order parameter, the step density, and a typical time scale, r s . 
The step density (re-)equilibrates locally at a time scale r s . The RSOS model is 
desribed by the KPZ equation only at time resolutions less than t 8 . In section 5, we 
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combine these results into a description of the global structure of the RSOS model 
phase diagram. We describe also how skewness can be introduced and controlled 
within the Langevin equation, eq.(l.l). 

Stationary skewness raises an important fundamental level-2 type issue. One 
of the basic premises of renormalization theory that all short-distance complexities 
are expressible in terms of irrelevant operators. Stationary state skewness should be 
an example of this. It is important to demonstrate this explicitly, since level-2 type 
understanding of KPZ type growth is still rudimentary. According to the empirical 
scaling theory, the growing interface displays scale invariance at large length scales 
and long times. It is invariant under the transformation: x' = bx, h' = b a h, and 
t' = b z t. The moments of the height distribution scale as 

W n (L, t) = b^Wnib^L, b~ z t) (1.3) 

a is the stationary state roughness exponent and z the dynamic exponent. Galilean 
invariance in the Burgers equation representation of eq.(l.l) yields the identity 
a + z = 2 [1-5,10]. In ID the height variable scales with the exponent a = 1/2. 
The width of the surface, W2, diverges linearly with L in the stationary state. The 
step-step correlation function 

G 8 (r) = {Vh x+r Vh x ) ~ exp(-r/60 (1-4) 

has a finite correlation length, or at least decays fast enough, such that going up 
and down along the surface becomes a random walk at large length scales. It is 
unlikely that KPZ type models with stationary state skewness belong to a different 
universality as those without it. The numerical evidence for z = 1.5 is very strong, 
in particular at the KK point. But, are they all described by the same fixed point? 

The conventional KPZ fixed point lies at zero skewness. It applies to eq.(l.l), to 
the BCSOS model, and also to the RSOS model inside the non-skewed 4 dimensional 
subspace. This fixed point describes skewed surfaces as well, if we can demonstrate 
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that the skewness crossover operator is irrelevant in the sense of renormalization 
theory. This does not mean that aspects like skewness vanish in the thermodynamic 
limit. Skewness W3 is allowed to diverge with system size, but at a rate slower than 
its naive exponent 3a. In section 5 we show numerically, using a master equation 
finite size scaling analysis, that the skewness in the RSOS model diverges in the 
thermodynamic limit, as W$ ~ N x , with x ~ 0.5. This confirms that the KPZ 
fixed point lies at zero skewness. The skewness crossover scaling exponent is equal 
to y sk ~ -1. 
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2. THE STATIONARY STATE IN THE BCSOS MODEL 



The stationary state of the BCSOS growth model is simple. It the completely 
disordered state. This is a consequence of the symmetry properties of the master 
equation. First we review these symmetries and then rederive the stationary state 
by a more complex method, one that can be generalized to the RSOS model in 
section 3. 

Consider a surface built from rectangular shaped bricks (a conventional brick 
wall turned over 90°). The surface heights at x = 2n + 1/2 are even integers and 
those at x = 2n — 1/2 are odd. Nearest neighbour columns differ in height by one 
unit. Each bond x = n contains a step S n = ±1. The growth rule is as follows: 
Choose one of the columns at random. If this column is at the bottom of a local 
valley, «S n _i = —1 and S n = +1, a particle adsorbs with probability p (and nothing 
happens with probability 1 — p). If it is at the top of a local hill, S n -i = +1 
and S n = —1, a particle desorbs with probability q (and nothing happens with 
probability 1 — q). If it is part of a local slope, S n -± = S n , nothing happens. 
This model has been studied extensively in the literature. First by Monte Carlo 
simulations [28,3]. More recently it was realized it can be solved exactly [6-9]. 

The Master equation 




(2.1) 



describes the time evolution of the probability distribution c({S n }) 




(2.2) 



{^} 



The time evolution operator has the familiar form 




(2.3) 



n 



but 7i is not Hermitian. 



n(n, n + l) = \e[l - S z n S z n+1 - 2(S+S~ +1 + S~S+ +1 ) 




(2.4) 
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with s = (p — q)/(p + q), and periodic boundary conditions, S^ +N = S^. T must be 
applied N times to evolve the system by one unit of time. Without loss of generality 
we can set e = p + q = 1. 

At s = 0, Ti is identical to the Hamiltonian of the so-called XXZ quantum 
spin- 1 chain. For any value of s the model is equivalent to the 2D equilibrium 6- 
vertex model in an electric field s. To be more precise, the time evolution operator 
is identical to the transfer matrix of the 6-vertex model when the sites are being 
updated sequentially instead of at random. That master equation reduces to eq.(2.3) 
in the time continuum limit. The 6-vertex model and eq.(2.4) are exactly soluble 
[7,8,31-33]. The master equation follows a special line through the 6-vertex model 
phase diagram, where T is a stochastic matrix (T preserves probability). The 6- 
vertex model describes the temperature evolution of equilibrium crystal surfaces 
[32]. KPZ type growth maps onto facet-ridge endpoints [13, 31], special points in 
the phase diagram where the rough (chiral-Luttinger liquid) and faceted phases 
meet. 

The stationary state of eq.(2.3) is very simple. It is the disordered state \D) 
where all coefficients c({S n }) are equal. For non-growing surfaces, with p = q, 
this is obvious. The dynamic rule is then equivalent to a Monte Carlo process 
in equilibrium statistical mechanics. The coefficients of the stationary state are 
proportional to the equilibrium Boltzmann weights. \D) is the equilibrium state 
since this ID BCSOS model lacks any interactions. 

\D) is the stationary state for the growing BCSOS surface as well. Stochastic 
processes preserve probability. Algebraically, this is expressed by the property that 
(D\ is the left eigenvector of the largest eigenvalue A = 1 for all stochastic time 
evolution operators. In general, the adjoint of a time evolution operator is not 
stochastic, but for the BCSOS model it is. describes the same growth process, 
but with the role of particles and holes reversed. The left and right eigenvectors 
switch role, and therefore \D) is the right eigenvector for all s^O. The stationary 
state is completely disordered. 
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The fact that is stochastic in the entire phase diagram is an accident. It is 
also an accident that is identical to the particle-hole transformed time evolution 
operator. We will see this clearly in section 3 during the discussion of the RSOS 
model. For that model we will be able to generalize the following method to derive 
the stationary state. 

The BCSOS model can be interpreted as a model for (classical) particles hop- 
ping along a chain in an electric field. S* = 1 represent an occupied site and 
S 1 ^ = -1 an empty site. This is also known as the asymmetric exclusion model 
[14]. Consider the probability c({S n }) to be in micro state {S n }. Let the dots in 
c(.., S* +l7 ..) mean that all spins other than those at sites n and n + 1 are the 
same for all states considered and do not change. The action of each individual 
H(n, n + 1), 

c(-.,+, -,-)t+i/N = (1 -q) c(.., +,-,- )t + P c(.. ,-,+,..)* (2-5) 
leaves \^>) invariant when 

+>->") = g (2.6) 
c(.., -,+,..) q 

The electrostatic energy of charged particles hopping along a wire in an electric 
field £ is equal to £/C, with 

n 

A conventional MC simulation rule for such a system obeys detailed balance, in 
which the transition probabilities between states a and b are related as 

= exp (5 AK/k B T). (2.8) 

This is the precisely the content of eq.(2.6). The stationary state is the "atmospheric 
law" density distribution: 



pin) = ex.p(—£n/kBT) 
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(2.9) 



and p/q = exp (£//c B T). This is not the solution we are looking for. It applies 
to the wrong boundary conditions. Eq.(2.9) is valid for open boundary conditions. 
In our case the particles keep running around in a circle and no density profile 
can build up to stop this flow. The detailed balance approach is apparently too 
restrictive. We must attack the problem less locally. 

The following property will prove essential. The number of hill tops and valley 
bottoms is the same in every configuration. This is easily established by drawing a 
typical configuration with periodic boundary conditions and a specific average slope. 
Algebraically it can be shown as follows. Consider the four nearest neighbour step- 
step densities in each configuration: , d |_, d ++ , and d They represent 

respectively, the density of hill tops, valley bottoms, and up and down slopes. The 
density of up-spins, d + , is related by normalization as 



This implies that cL| = d |_. 

Let's test the following assumption: The probability distribution c({S n }) de- 
pends only on the total number of particles Af and the total electrostatic energy 
/C. The average slope of the surface is fixed by the boundary conditions. Therefore 
Af is a constant of motion and can be ignored. This leaves only the dependence on 
the electrostatic energy, c({S n }) = c(/C). This probability distribution time evolves 
according to eq.(2.3) as 



The first term on the right hand side represents the contributions for when nothing 
happens between t — > t + 1/N. The second and third term represent events where 
one particle adsorbs or evaporates. Adsorption destroys one valley bottom and 
creates one hill top. Evaporation does the opposite. The adsorption probability 
is proportional to <i_| since this could have happened to any of the hill tops in 



d + = d ++ + d- + = d ++ + d + - 



(2.10) 



<9c(/C) 
dt 



(pd |_ + qd-\ )c(X) + pd^ c(/C + 1) + qd |_c(/C — 1) 



(2.11) 
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existence at time t + 1/N. It is also proportional to C(/C + 1), because in the 
electric field interpretation one particle hops in the direction of the electric field. 

Our assumption for the stationary state is correct if the right hand side of 
eq.(2.11) vanishes for all dij. This seems to imply two conditions, but there is only 

one, since d |_ can be eliminated using the identity d |_ = d-\ The right hand 

side vanishes for all configurations if 

p + q=p/w + qw (2-12) 

with w = c(/C)/c(/C + 1). This yields two stationary state solutions: the detailed 
balance solution with w = p/q, and the stationary growing state solution \D) with 
w = 1, the state we are looking for. 

Without the global identity d |_ = d^ , eq.(2.10) would yield two conditions, 

one for d |_ and one for cL| These reproduce only the detailed balance solution. 

The stationary growth solution with w = 1 is the result of the non-local property 
that there are just as many valley bottoms as hill tops. Every hill top (valley 
bottom) could have been created during t — > t + 1/N with probability p (q) and 
every valley bottom (hill top) can be destroyed with the same probability. Since 
d- + = d + - the right hand side of eq.(2.11) vanishes when the coefficients c{K) are 
independent of /C. 
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3. THE STATIONARY STATE IN THE RSOS MODEL 



The RSOS growth model describes the growth of simple cubic surfaces in which 
only mono-atomic steps are allowed. Nearest neighbour columns can differ by only 
dh = S* = 0, ±1. The master equation is more complex than eq.(2.3) for the 
BCSOS model, since more types of local configurations are possible. Each has 
its own transition probabilities. After choosing one of the columns at random, one 
particle can be evaporated or deposited at x = n+ \ with a probability that depends 
on the height differences and S* +1 . The time evolution operator is again of the 
form eq.(2.3) with 

n= £ {[ Ph (i-s+s~ +1 ) + qv (i-s-s+ +1 ) ] s(o) n 5(o) n+1 

n=l,N 

+ q h (1-S~S+ +1 ) 8(+) n S(-) n+1 

+ Pv (l- S+S~ +1 ) 6(-) n 6(+) n+1 

+ Ps (1-S+S~ +1 ) [S(0) n 5(+) n+1 +5(-U(0) n+1 ] 

+ q s (l- S~S+ +1 ) [8(+) n 8(0) n+1 + S(0) n 5(-) n+1 ] } (3.1) 

and the following definitions: 5(0) = (1 + S z )(l - S z ), 5(+) = \S Z (S Z + 1), and 
5(—) = ^S Z (S Z — 1). The raising and lowering operators S + and S~ are normal- 
ized such that S z = S+S- - S~S+ and S+ S~ + S~ S+ = 2 - (S z f . This time 
evolution operator has 6 parameters. They have the following interpretations: ph is 
the adsorption probability and q v the evaporation probability when the surface is 
locally flat; qt is the evaporation probability at a local hill top; p v is the adsorption 
probability into a local valley bottom; p s is the adsorption probability and q s the 
evaporation probability at a step. 

7i is not Hermitian, but resembles the Hamiltonian of a quantum spin-1 chain. 

H D (n,n + 1) = c s + (sm - 2c s )(^) 2 - \a s S z n S z n+l + (c s - s m + \a s ){S z n S z n+1 ) 2 
+(h s d — jf l g)S z S z +1 (S z +1 — S z ) (3.2a) 
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the diagonal part, and 



H OD (n, n+l)= -J^+S+S- +1 - F_S~ S+ +1 (3.26) 

the off-diagonal part, with 

T + =(PH- Ps )(S z n S z n+1 f + (p v - Ps )[l - (5*) 2 ][1 - (S z n+1 ) 2 ] + Ps 

T~ =(Qv ~ q s )(S z n S z n+l ) 2 + (q h - q s )[l - (5*) 2 ][1 - (^ +1 ) 2 ] + q a . (3.2c) 

In Hd we introduced a second notation for the 6 parameters: c s = pn+Qv is the step 
creation probability; a s = p v + qh is the step annihilation probability; s m = p s + q s 
is the step mobility; f g =Ph~ Qv is the growth probability at flat surface areas; 
h g = p v — qh is the growth probability at hill tops and valleys; and Sd = p s — q s is 
the growth probability at steps (step drift) . 

The phase diagram is only 5 dimensional because rescaling all 6 parameters by a 
common factor redefines the unit of time. Some familiar models are contained in this 
master equation as special points or lines [29,30]. For example, the Kim-Kosterlitz 
model [29] corresponds to q v = qh = q s = (no evaporation) and p v = Ph = Ps- 

Recall (see section 2) that the stationary state of the BCSOS model, eq.(2.4), 
is trivial. It is the completely disordered state \D), because the adjoint of its time 
evolution operator it also stochastic. It describes the time-evolution of the same 
surface in the particle-hole transformed representation. The adjoint of eq.(3.1) is 
only stochastic in a 3D subspace of the phase diagram. The diagonal parts of T' 
and T are identical, but the coefficients in the off-diagonal part switch position, 

7~(c s , <x s , s m ; fg, h g , Sd}^ = 

Cs-i &si S m 

; fg, hg, s d ) + T OD (a s , c s , s m ; -h g , -f g , -s d ) (3.3) 

The parameter space in which preserves probability is larger than the 2D self- 
adjoint subspace in which the surface does not grow (ph = qh, Pv = 1v, and p s = q s ). 
7~t is stochastic whenever the diagonal part To in eq.(3.3), is invariant under the 
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transformation c s <-> a s h g <-> — / s , and — > — (the way the parameters switch 
in the off-diagonal part). This is true for 

c s = a s , h g - 2s d = -f g + 2s d . (3.4) 

Inside this 3D subspace the adjoint of the time evolution operator is stochastic and 
describes surfaces in which the probabilities are interchanged as 

T"(C S , <X S , S m ', fg, hg, Sd)^ = T~{tt.Sl ^771 ? hg, fg, S (f) (^-5) 

In the BCSOS model, 7~t is identical to the particle-hole transformed dynamics. 
This is not true anymore in the RSOS model; the particle-hole transformed time 
evolution operator 

Oph T{c s , a s , s m ; f g , h g , s d ) = T(c s , a s , s m ; -f g , -h g , -s d ) (3.6) 

is different from the one in eq.(3.5). 

Next, we generalize this stationary state to a 4D subspace. The BCSOS model 
master equation describes a lattice gas of non-interacting charged particles in an 
electric field. The RSOS master equation is a particle-hole generalization of this. 
The Si^ = ±1 states represent sites occupied by particles with electric charge ±1. 
S* = states represent empty sites. The energy of such a lattice gas is of the form 

E = Yy»{S z n ? + £n S* n \. (3.7) 

n 

The chemical potential \i is needed, because the dynamic rule does not conserve 
the total number of particles (the number of steps in the surface). The chemical 
potential for the total charge can be omitted because electric charge is conserved. 
We are looking for a stationary state of the form c({S^}) = c(A/" s ,/C) in which the 
coefficients only depend on the number of particles, Af s and the total electrostatic 
energy, K. = ^n^, 

<{Sn}) ~ exp[-//jV s - £K] = z^ Ns w K . (3.8) 
15 



Let us try the conventional detailed balance approach for Monte Carlo simulations. 
Require that the action of each individual H nj7l+ i leave \^) invariant. This yields 
the conditions 

q h c(..,0,0,..) 



Ph c(..,+, -,..) 
p v c(..,0,0, ..) 



Qv c(.., -,+,..) 

Ps_ _ c(..,+,0, ■■) = c(..,0, -, ..) 

q s c(..,0,+, ..) c(..,-,0, ..) 



(3.9) 



(using the same type of notation as in eq.(2.5)). They are satisfied by the stationary 
state of the form eq.(3.8) when 

Ph/Qh = exp[(£ - 2n)/k B T] 
Pv/Qv = exp[(£ + 2p)/k B T\ 

p s /q s = exp[£//e B T] (3.10) 

This defines a 4D subspace inside the 5D phase diagram, since each set of values S 
and leaves two ratios unspecified, pn/Pv and Ph/p s - 

Eq.(3.9) is not the stationary state we are interested in. It applies to the wrong 
type of boundary condition. Let's apply the same non-local approach as in section 
2. Assume that the stationary state is of the form eq.(3.8) and express the time 

evolution in terms of the 9 nearest neighbour step-step densities <ioo> <^h )•••) d++. 

The equations of motion read: 

d C( ^'^ = - [{ph + q v )d o + p v d-+ + qhd+-] c(J\f s , K) 

+ Ph d+- c(Af s - 2, K + 1) + q v d-+ c(Af s - 2, K - 1) 

+ q h d 00 c{N s + 2, K, - 1) + p v d 00 c{N s + 2, K + 1) 

+ q s (d 0+ +d- ) c(K,IC-l)+p s (d + o + do-) c(A/" s ,/C + l) 

- [p s (d + + d- ) + q s (d +0 + d -)] c(Af s , /C) (3.11) 
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This is the generalization of eq.(2.11). The stationary state is of the form eq.(3.8), 
if the right hand side vanishes for all dij. Naively this yields 5 equations with 5 
unknowns. However, the dij in eq.(3.11) are not independent. The analogues of 
eq.(2.10) are 

d + =d ++ + d +0 + d + - = d ++ + d + + d-+ 

d_ =d +d -+d + -=d + d- + d- + (3.12) 

with d + and d_ the density of up- and down-steps, and do + d + + d_ = 1. Adding 
and subtracting these two equations yields two identities: 

2d + - + d +0 + d - =2d- + + d 0+ + d- 

d+o - d - =d 0+ - d-o (3.13) 

After employing the first one, only doo ? d^ , d |_, and d+o+do- remain in eq.(3.11). 

Eq.(3.8) is the correct stationary state if the following 4 conditions are satisfied: 

(Ph + Qv) =qhzw~ 1 +p v zw 

Qh = PhZ~ 1 w -2(p a - qsiv' 1 ) 
Pv = qvizw)' 1 + 2(p s - q s w~ x ) 
= ip s - QsW' 1 )^ - w) (3.14) 

The last equation indicates that there are two solutions: one with q s /p s = w and 
one with w = 1. The first solution, with w = q s /p s , reproduces the atmospheric 
law type charge profile (3.9). The second solution, with w = 1, applies to periodic 
boundary conditions and is the one we are looking for. For w = 1, the last condition 
in (3.14) is satisfied without any constraints on p s and q s . Only 2 of the 3 remaining 
conditions are independent. One sets the step density and the other specifies 1 
condition between the 5 growth parameters. The stationary state is of the form 

c({S z n }) = c{Ns) ~ exp[-^AT s ] = (v 7 ^ (3.15a) 
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inside the 4 dimensional subspace denned as: 

Ph + Qv 



qh +p v 

= z(p v - q h ) + (p h - q v ) - 4z(p s - q s ) (3.156) 
the latter is equivalent to 

z = c s /a s 

4s d z = zh g + f g (3.15c) 

Both solutions of w are valid in a subspace of co-dimension 1 of the 5 dimensional 
parameter space, but these subspaces do not coincide, nor do the equations for the 
step-fugacity z. 

In the stationary state (3.15) the step density is equal to 

^ = ms z n ?m = y^j= z (3.i6) 

for zero surface tilt, v = (D\S*\S) = 0. Notice that for z = 1 we recover the 
disordered state solution, eq.(3.4). 

The stationary growth rate, r g , is equal to 

r g = (D\(J^ + S+S- +1 -^S-S+ +i m 

=f g (l-s) 2 + ±h g s 2 + s d s(l-s) (3.17) 

which simplifies with eq.(3.15c) and eq.(3.16) into 

r g = s d s (3.18) 

The height-height correlation function diverges linearly because the spins are un- 
corrected: 

(D\(h n+h - h n+mH r\y) = (d\ Yl (^) 2 i*) = ms - ( 3 - 19 ) 

n<i<n+m 
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This confirms that the surface roughness critical exponent is equal to a = ■§. 

It is amazing that the stationary state is this simple in such a large fraction of 
the phase diagram. The step-step correlation function, eq.(1.4), has a sharp cut-off, 
to such an extent that steps are completely disordered. Outside this 4D subspace 
the stationary state becomes more complex. Eq.(1.4) does not have a sharp cutoff 
anymore. We checked this numerically. Attempts to extend the above derivation to 
a form in which the c({S n }) are functions only of the nearest neighbour step-step 
correlations are doomed. The number of equations increases rapidly and they can 
not be satisfied inside the 5D parameter space. 
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4. MEAN FIELD THEORY AND LANGEVIN EQUATIONS 



In this section we apply mean field theory to the BCSOS and RSOS model. The 
mean field equations of motion for the order parameters reproduce the deterministic 
part of the KPZ equation. This is similar to equilibrium critical phenomena, where 
mean field theory for the Ising model reproduces only the non-fluctuating part of 
4 -theory. The object is to better understand the structure of the 5D RSOS model 
phase diagram; in particular, to relate qualitatively the parameters of the RSOS 
model to those in the KPZ equation. 

Consider the BCSOS Master equation, eq.(2.4). The equation of motion 



dv 
dt 



- = v n (t + 1) - v n (t) = (D\[SZ,T\\*) t (4-1) 



for the local slope of the surface 

vn(t) = (D\SZ\9) t (4.2) 

does not close because of the presence of nearest neighbour step-step correlations 
on the right hand side. Define 

d(a, b) x = {D\5{S z n - a)5(S* n+1 - b)\9) (4.3) 

with a, b = ±1 and x = n + |, and replace the slope variable by the surface height 
using the relation v n = h x+ i — h x . Eq.(4.1) yields, 

— =pd(+,-) x -qd(-,+) x (4.4) 

This is the top of a BBGKY hierarchy of equations of motion. In lowest order 
approximation all step-step correlation are neglected. The d(a, b) x are assumed to 
factorize. The equation of motion for the average local slope v n (t) then closes. 

d 1 1 

-Q t K = -p{l - v n )(l + v n+1 ) - -q(l + v n )(l - v n+1 ). (4.5) 

20 



A faster route to exactly the same equation is to apply the mean field approximation 
directly to the master equation. The adsorption and evaporation transition rates 
at site x are assumed to depend on the expectation values of the step density v n 
instead of their actual values S*. Eq.(4.5) is identical to the deterministic part of 
the KPZ equation, eq.(l.l), in the spatial continuum limit 

dh 1. , 1. ,d 2 h 1. w<9/i,o „ N 

at = 4<p-«) + i<p+»>flii-i(i>-»>(flS? (46) 

A is negative and reduces the growth rate at slopes. This reflects that in the BCSOS 
model sloped surface sections, with S n = S n +i, are inactive. 

Apply the same type of approximation to the RSOS model. The first observa- 
tion is that the RSOS model contains two point-like expectation values instead of 
one: the local step density, s n , and the local slope of the surface, v n , 

Sn(t) = (D\(S* n ) 2 \*) t 

vn(t) = (D\S z n \*) t . (4.7) 

Mean field theory involves two order parameters instead of one. The equations of 
motion for these two order parameters, 

- = 8i (t + 1) - Si (t) = (D\[(S?f, T\\*)t (4.8.a) 

^ = Vi (t + 1) - Vi (t) = (D\Sf, m*> t (4.8.6) 

are again the top of a BBGKY type hierarchy of coupled differential equations: 

ds r , i 

=(Ph + q v )[d(0, 0) x + d(0, 0) x+1 ] 

-p v [d(-, +) x + d(-, +) x+ i] - qh[d{+, -) x + d(+, -) x +i] 

+ Ps [d(-, 0) x - d(0, +) x - d(-, 0) x+1 + d(0, +) x+1 ] 

+ q s [d(+, 0) x - d(0, -) x - d(+, 0) x+1 + d(0, -) x+1 ] (4.9.a) 



dt 



= (Ph ~ q v )d{0, 0) x + p v d(-, +) x - q h d(+, -) x 

+ Ps [d(-, 0) x + d(0, +) x ] - q s [d(+, 0) x + d(0, -) x ] 
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(4.9.6) 



In mean field theory the nearest neighbour step-step correlations factorize. The 
resulting equations are not very transparent. Therefore, we first consider the special 
case where the surface is spatially uniform, such that s = s n and h = h x do not 
depend on position. The mean field approximation for eq.(4.9) then reads 
ds 

— =2c s (l - s) 2 - ±a s (s - v)(s + v) (4.10.a) 

-gf =fg(l - S) 2 + \h g {s 2 - V 2 ) + 8d (l ~ 8)8 (4.10.6) 

The equation for the step density contains a "mass" term, which defines a char- 
acteristic time T a . At zero net tilt, v = 0, the step density reaches its stationary 
value 

s = [l + iJ-)]" 1 ( 4 - 1La ) 



V c s 

after a characteristic time 

r^ 1 = 2v^7 (4.11.6) 
At time scales larger than r s the surface grows at a constant average rate 

r = f g (l - s ) 2 + \h g sl + s d (l - s )s (4.11.c) 

This means that, although there are two order parameters, only one of them, the 
local slope v n , fluctuates at time scales larger than r s . The full mean field equation 
for the step density reads: 

^ = 2c s (l - s) 2 - ±a s (s 2 - v 2 ) + s d |-[(l - s)v] + hh 9 (v^ - s^) 



,d 2 s d 2 v . . d 2 s 



in which all derivatives are discrete. Define A as the deviation from the stationary 
value, s = s + A. At time scales larger r s the step density does not behave as 
an independent dynamic variable. A follows local fluctuations in the slope of the 
surface such that the right hand side of eq.(4.12) remains equal to zero: 
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with 5 = ±a s and e = s<i(l — sq) ~~ h^g s o- The full equation of motion for /i x 
is complex. It contains many terms, involving the step density, derivatives of the 
surface slope, and combinations of these and their derivatives. At time scales larger 
than t s , the step density is not an independent variable, and can be eliminated 
using eq.(4.13). The equation for dh/dt then reduces to the deterministic part of 
the KPZ equation, eq.(l.l), with 

A ~ ±a s T s [h g s - 4/ s (l - s ) + 2s d (l - 2s )] - \h g , 

v ~ \s m {\ - s ) + ia s s . (4.14) 

and many higher order terms as well. The above two derivations illustrate that 
the BCSOS and RSOS growth models belong to the KPZ universality class. That 
is hardly a surprise however. The significance is the identication of an additional 
order parameter and of the time scale r s . At short time scales the step density 
behaves as an independent dynamic variable and the KPZ equation description is 
incomplete. Fortunately r s is typically very short, of the order of only a few time 
steps. 
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5. SKEWNESS IN THE RSOS MODEL 



The phase diagram of the RSOS model seems too large to visualize, but the 
results of the previous sections put some order into it. The uncorrelated stationary 
state (in the subspace defined by eq.(3.15)) and the mean field theory define a set 
of characteristic densities and time scales. The six parameters in eq.(3.2) fall into 
two groups: c s , a s , and s m , are equilibration type parameters; f g h g and Sd, are 
growth type parameters. 

The ratio between the step creation and annihilation probabilities, z = c s /a s , 
controls the steps density, s. Their product r s _1 = 2^a s c s controls the time scale 
at which the step density equilibrates. The step mobility s m is a suitable unit of 
time. (We set s m = 1.) KPZ type scaling is realized at time scales larger than r s 
and length scales larger than the inverse of the step density. This leaves us with an 
effective 3D phase diagram, characterized by the growth parameters, f g , h g , and 
Sd- The linear combination 

r g = f g (l - s 2 ) + \h g s 2 + s d s(l - s), (5.1) 

with s defined in eq.(3.15), controls the growth rate. The linear combination 

u s k = zh g + f g - Azs d (5.2) 

controls the skewness of the stationary state, see eq.(3.15). In the zero skewness 
plane, u s k = 0, r g is exactly equal to the growth rate and s is exactly equal to 
the step density. The third independent combination of f g h g and Sd controls the 
strength of the non-linear term of the KPZ equation. Eq.(4.14) is an approximation 
for A. 

The phase diagram of the RSOS model is effectively only three dimensional, but 
is still larger than the one for the KPZ equation, eq.(l.l). The RSOS model allows us 
to control skewness. Stationary skewness can turned on in the KPZ equation as well. 
The stationary state of eq.(l.l) is known exactly in ID. It is a Gaussian distribution 
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without any skewness [5,10,27]. This is an accident. Additional operators in the 
equation should introduce stationary state skewness. The following A2 type term is 
an example of this 

^=r 9 + uV 2 h + ±A(V/i) 2 + A 2 (V 2 / i ) 2 + V . (5.3) 

We performed some (qualitative) Monte Carlo runs at A 2 7^ [34]. Those confirm 
that the stationary state is skewed. The BCSOS growth model and the u s k = 
RSOS model represent special cuts through the extended KPZ equation, where 
parameters such as A2 take special values such that the stationary state lacks skew- 
ness. On a qualitative level the role of A2 can be understood as follows: In a flat 
surface, the growth rate in the RSOS model is determined by f g , and in the KPZ 
equation by r g . In a sloped surface, the growth rate in the RSOS model is deter- 
mined by f g (l — so) + SdSo, and in the KPZ equation by r g + A. In a hilly surface 
without skewness with densely packed steps, the growth rate in the RSOS model is 
determined by h g , and in the KPZ equation, eq.(5.3), by r g + A 2 . 

We hoped that mean field theory (section 4) would provide a meaningful es- 
timate of the u s k = space in the generalized KPZ equation. Unfortunately, too 
many operators seem to be involved. For example, the analysis of section 4 gives a 
non-zero value for A 2 in the BCSOS model. Apparently that value for A 2 is com- 
pensated by significant contributions of other operators such that stationary state 
skewness remains absent. 

The reverse route might be more promising. The major difference between 
eq.(5.3) and microscopic models like the BCSOS and RSOS model is the manner 
in which the noise couples to the local surface structure. In the Langevin equation 
they are "additive". The probablity distribution obeys a Fokker-Planck equation 
[5,35]. For example, in discretized time it is allowed to visualize the time evolution 
as determinsitic between t — > t + \ and purely stochastic (ballistic deposition like) 
between t + \ — > t + 1. 

h(t + i) =h(t) +r g + uV 2 h + iA(V/i) 2 + A 2 (V 2 /i) 2 
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h(t+l) =h{t+±) + rj. 



(5.4) 



In the BCSOS and RSOS model the noise couples much more directly to the local 
surface configuration. Correlation effects between steps are more pronounced. 

The KPZ fixed point must lie inside the non-skewness subspace. As a rule, 
fixed points lie in subspaces where the short distance properties mimic as much as 
possible the proper long distance scaling properties. For example, the fixed point 
of the Ising model on a square lattice lies in the subspace where the next nearest 
neighbour interactions are of the same order of magnitude as the nearest neighbour 
interactions, such that the correlation functions are rotational invariant not only at 
large distances, but also at short distances. The cubic anisotropy is an irrelevant 
operator and needs to be turned-off at the fixed point. 

The following numerical results demonstrate that stationary state skewness is 
indeed irrelevant in KPZ type growth. We determine the largest eigenvalues and 
eigenvectors of the time evolution operator for system sizes 4 < L < 14. These val- 
ues are exact, and allow a detailed finite size scaling (FSS) analysis. This method is 
identical to conventional transfer matrix FSS calculations in equilibrium phase tran- 
sitions, and works equally well for master equations [13]. We applied this method 
to several points in the phase diagram, but we present only our results for the Kim- 
Kosterlitz (KK) model [29]. This is a typical point outside the u s k = subspace. 

On a local level skewness implies an imbalance between sharp hill tops and 
valleys bottoms. The quantity 

Psk ~ d(+,-)-d(-,+) (5 - 5) 

measures the density difference between sharp hill-tops and sharp valley-bottoms. 
p s k converges to a non-zero value at the KK point, as shown in Fig. 2. This illustrates 
the presence of skewness in the stationary state at short distances. The long distance 
probe for skewness is the third moment, W3, of the height distribution function, see 
eq.(1.2). 
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The scaling behaviour of W 3 tells us whether stationary state skewness is 
present or absent at the KPZ fixed point. It should diverge as W3 ~ L 3a if present. 
Otherwise, we must incorporate into eq.(1.3) the crossover scaling in the skewness 
direction u sk , 

W n (L,t,u sk )=b na W n (b- 1 L,b- z t,b y ° k u sk ). (5.6) 
and expand this in small u sk , 

W n (L, 00, u sk ) ~ L na [W(l, 00, 0) + LV*u ak W'(l, 00, 0) + ....] (5.7) 

In the absence of skewness the amplitude W^l, 00, 0) of the leading term is equal 
to zero and the third moment scales as W3 ~ AL X with x = 3a + y sk . Fig. 3 
demonstrates that stationary state skewness is indeed absent at the KPZ fixed 
point. The exponent x = 0.4 ± 0.1 suggests that the skewness crossover exponent 
is equal to y sk = —1 (the nearest integer). The amplitude A of the third moment, 
see Fig. 4, must be proportional to u sk . The local measure for the skewness, p sk , 
in Fig. 2 should be approximately proportional to u sk as well. We calculated the 
ratio A/p sk at several points in the phase diagram, and find it indeed to be almost 
a constant. 

In conclusion, in this paper we study the phase diagram of the RSOS growth 
model in 1+1 dimensions. Its phase diagram contains a 4D subspace in which the 
stationary state is completely uncorrelated. Familiar models, such as the Kim- 
Kosterlitz model, lie outside this subspace. Their stationary states contain addi- 
tional features, like skewness. Stationary state skewness diverges with system size, 
but such that the scaling properties are still described by the KPZ fixed point at 
zero skewness. The skewness crossover exponent is equal to y sk ~ —1. 

This research is supported by NSF grant DMR-9205125. 
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FIGURE CAPTIONS 



1. Temporal skewness in the deterministic KPZ equation. Time evolution of a 
sinusoidal initial state. 

2. Finite size scaling of the local skewness order parameter, p s k, eq.(5.5), at the 
Kim-Kosterlitz point. p s f- represents the density difference between sharp hill 
tops and sharp valley bottoms. 

3. The scaling exponent x of the third moment, W3 ~ AN X , at the Kim-Kosterlitz 
point. The drawn line represents finite size scaling estimates x{N) from the 
exact values of W 3 (N) at N-l and N+l for N = 5, 7, ..13. The dashed line 
shows estimates for x defined as x(N) ~ x + B/N 2 at successive values of N. 

4. Finite size scaling estimates of the amplitude A of the third moment, W% ~ 
AN X , at the Kim-Kosterlitz point, assuming that x = 0.4. The dashed line 
represents the same type of extrapolation as in Fig. 3. 
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